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We present a structural analysis of the graphene-4HSiC(0001) interface using sur- 
face x-ray reflectivity. We find that the interface is composed of an extended recon- 
struction of two SiC bilayers. The interface directly below the first graphene sheet 
is an extended layer that is more than twice the thickness of a bulk SiC bilayer 
(~ 1.7A compared to 0.63A). The distance from this interface layer to the first 
graphene sheet is much smaller than the graphite interlayer spacing but larger than 
the same distance measured for graphene grown on the (0001) surface, as predicted 
previously by ab intio calculations. 

PACS numbers: 68.55.-a, 68.35.-p, 61.10.Kw, 61.46.-w 
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I. INTRODUCTION 



The direct growth of single or multiple sheets of graphene on an insulating or semi- 
conducting substrate is known as epitaxial graphene (EG). Because this material has been 
identified as a viable all-carbon candidate for post CMOS electronics, [l| there is a strong 
impetus for the study of EG and how it can be produced, lithographically patterned and 
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made into electronic devices. [2|, y|, |4J, [5|, la] The current substrate for EG growth is either 
of the two polar faces of hexagonal SiC: the SiC(OOOl) Si-terminated surface (Si-face) and 
the SiC(OOOl) C-terminated surface (C-face). Multilayer epitaxial graphene films grown on 
C-face SiC show electronic properties expected for an isolated graphene sheet including a 
Berry's phase of 7r, weak anti-localization, and a square root dependence of the Landau level 
energies with applied magnetic field. 0, 4, 3, Is,!^] The fact that transport in graphene grown 
on SiC substrates is so similar to transport properties expected for an isolated graphene sheet 
is quite remarkable since graphene/substrate interactions and multilayer graphene stacking 
should potentially influence the 2D Dirac electrons responsible for the unusual properties of 
graphene. On the other hand, these findings are very fortuitous and have made graphene 
grown on SiC the focus of research targeting a path towards graphene electronics. 

Studies of the graphene/SiC substrate and multilayer graphene stacking have begun to 
elucidate the relationship between EG's electronic properties and its structure. It is very 
clear from the earliest measurements that both the growth and structure of C-face and 
Si-face grown graphene are very different. 0, [h], 
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Graphene prepared on SiC 
in ultra high vacuum (UHV) grows relatively slowly on the Si-face as compared with the 
11| and tends to be ~ 1 — 5 layers thick, while graphene grown on the C-face grows 



C-face 



rapidly at lower temperatures compared to the Si-face and can achieve thickness well above 

n n n 

5 layers. [5|, [12|, [13( In addition, graphene grown on the C-face in UHV contains a high 

16| while graphene growth on 



concentrations of graphene nanocaps 
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and other defects. 



the Si-face causes the SiC substrate to roughen. [5| Furnace grown C-face graphene on the 



other hand reaches domain sizes > 5/im, 17|] many orders of magnitude larger than those 
grown in UHV on either the C-face or Si-face. [3] Correspondingly, electron mobilities and 
electronic coherence lengths are an order of magnitude higher on the furnace grown C- 
face graphene compared to UHV grown Si-face graphene. 3) J|j] Finally, Si-face graphene is 
epitaxial with (6\/3 x 6v / 3)R30 periodicity as observed by low energy electron diffraction 



(LEED) (i.e, the graphene is rotated 30° relative to the SiC < 1010 > direction) 
C-face films, on the other hand, can have multiple orientational phases 1(3, 
a unique stacking order. 
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Because of the relative ease of growing graphene in UHV on SiCfOOOl), the Si 
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y studied of the two hexagonal polar surfaces. |19|, |20J, 
, 32| Despite this body of work, the structure of this interface is still 
unknown even though it is crucial for understanding the electronic structure of subsequent 
graphene layers and the charge transfer between the substrate and the graphene film. Early 
studies of the graphitization of the Si-face presumed that the (6^x6^)1^30 pattern was due 
to the commensurate alignment of a graphene overlayer with an unreconstructed SiC surface, 



forming a moire pattern. [IOJ, [20|, [2l|, |22j However, the consensus from recent experiments by 
many groups is that the (6^x6-^/3)1^.30 structure is a true structural precursor phase to 
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brmation that persists, although slightly altered, after graphene has formed. 23j, 



33| While scanning tunneling microscopy (STM). 
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outlined some of the lateral features of the (6a/3 x 6v / 3)R30 structure, the buried interface 
structure remains elusive. 

In this paper, we use specular x-ray reflectivity to measure the buried structure of the 
graphene/4H-SiC(0001) interface. We find that the first layer of carbon with an areal density 
of graphene sits close to the last layer of atoms in the interface layer. For the Si-face the 
graphene-interface spacing is found to be 2.32 ± 0.08A. This number is consistent with 
ab-initio calculations that predict a covalently bonded insulating first graphene layer that 
"buffers" subsequent graphene layers from the substrate. 36|, |37J However, unlike the simple 
bulk terminated interface used in these calculations, the interface layer is found to be more 
complex. The structure of this interface layer suggests it plays an important role mediating 
the interaction of the SiC substrate with the graphene film. 



II. EXPERIMENTAL 



The substrates used in these studies were 4H-SiC purchased from Cree, Inc. 
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Prior 



to graphitization the 3 x 4 x 0.35mm samples were ultrasonically cleaned in acetone and 



ethanol. The Si-face samples were H 2 etched 171] and subsequently graphitized in UHV 
(P < 1 x 10~ 10 Torr) by electron-bombardment heating. Substrates were first heated to 
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1100°C for 6 min to form a (v3 x y / 3)R30 reconstruction, after which they were heated to 
1320°C for 8 min to remove surface contamination and to form a well ordered (Qy/3 x6-\/3)R30 
reconstruction. The samples were subsequently heated to 1400-1440°C for 6-12 min to create 
graphene films 1-3 layers thick. During graphitization the pressure in the system reached 
P < 1 x 10 -8 Torr. Despite the high pressure during growth, sample order is indistinguishable 
from Si-face sample order reported by other groups who grow at pressures ~ 10 _10 torr or 
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29j Growth-induced 



from those who use a Si flux to reduce surface oxides, 
substrate roughening in all reported investigations of Si-face graphitization leads to average 
SiC terrace sizes of ~ 300 — 500 A (graphene order is actually bigger than the terrace size 
because graphene grows over the SiC steps. [30j) It appears that the long range order of UHV 
grown graphene on the Si-face is not very sensitive to most details of the surface preparation. 

Once the samples were graphitized, they remained inert allowing them to be transported 
into the separate x-ray scattering chamber for analysis. The x-ray scattering experiments 
were performed at the Advanced Photon Source, Argonne National Laboratory, on the 
6IDC-yuCAT UHV (P < 2 x 10~ 10 torr) beam line. Experiments were performed at a photon 
energy of 16.2 keV. The number of graphene layers on the sample was determined by both 
Auger electron spectroscopy (AES) and x-ray reflectivity (as described in the next section). 
The reflectivity data is presented in standard SiC hexagonal reciprocal lattice units [r./.tz.] 
(h, k, £). These are defined by the momentum transfer vector q= (/iag iC , kb§ iC , £c§ iC ), where 
Og iC = 6g iC = 27r/(asicV / 3/2) and Cg iC = 27r/csic) [ fl SiC, &SiC and csic are the standard 4H-SiC 
hexagonal lattice constants] . The measured lattice constants were agjc = 3.079 ± O.OOlA, 



c SiC = 10.081 ± .002A and are within error bars of published values. 39J For reference the 



nominal hexagonal lattice constants for graphite are = 2.4589A, c G = 6.708A.|40| 



III. RESULTS 

To obtain detailed information about both the graphene films and the SiC-graphene 
interface, we have measured the surface x-ray specular reflectivity from graphitized 4H- 
SiC(0001). Details of the data collection and the model are similar to those used to determine 
the structure of the 4H-SiC(0001)/graphene system. [4l|] The data is collected by integrating 
rocking curves for different perpendicular momentum transfer vectors, q± = 2n£/csic, where 
q = kj — kj. Since the reflectivity only depends on q± (or equivalently £), the data can be 
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FIG. 1: A schematic model of multi-layer graphene grown on the 4H-SiC(0001) substrate. Dashed 
lines are the bulk SiC lattice planes before interface relaxation (A's). The 5 th plane of atoms 
(adatom) is displaced d a d from the topmost atom plane in the interface. (•) are carbon atoms and 
(o) are silicon atoms. The shaded circles in the interface ("layer-0") can be either carbon or silicon 
atoms. The graphene layers above the interface layer are referred to as "layer-", -2, -3, etc. 

analyzed using a one- dimensional model where all lateral information is averaged over the 
0.4 x 0.4mm x-ray beam. 

The schematic model of the graphene covered SiC(OOOl) surface is shown in Fig. [TJ This 
general model allows us to explore a number of possible graphene/SiC structures proposed 
by previous experiments. In the model the SiC substrate contribution is broken into two 
terms: (i) the amplitude from a bulk terminated surface, (ii) the combined amplitude from 
a relaxed SiC bilayer just above the bulk and a reconstructed interface layer (defined as the 
interface "layer-"). The scattered x-ray intensity 1(6, £) is then the result of a sum of three 
scattered amplitudes; the bulk Fb u ik, the interface region Fj, and the graphene Fq: 



1(6,. 



A(9,f)e- 47SiCSin2rf/2 



A(e, 



1 - e -»« psic 
is a term that contains all corrections due to the experimental geometry. 



(1) 
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43|, 



44J The exponential term accounts for the substrate roughness caused by half-cell step 

/2). 7sic 



fluctuations in the SiC surface (the predominant step height on 4H samples; [5] 
is the variance in the number of half-cell layers in the surface due to steps. 
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451 ] Roughly, 
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7sic is proportional to the SiC step density. The first term in Eq. (pQ) is the bulk 4H- 
SiC structure factor, i 7 buik(^),[46| modified by the crystal truncation term, (1 — e - 2me -}~ 1 
(Ref.|43). F G (£) in Eq. flU) is weighted by the ratio of the areal densities of a 4H-SiC(0001) 
and a graphene (0001) plane; pc/psi=3.132, to properly normalize the scattered amplitude 
from the graphene layer per 4H-SiC(0001) (1 x 1) unit cell. 

Fj(£) in Eq. ([1]) is the structure factor of the interface region between the bulk and 
the graphene film. Although we cannot obtain lateral information about the SiC(OOOl) 
(6v / 3x6\/3)R30 structure from reflectivity data, the vertical shifts of atoms and layer density 
changes associated with them can be determined. To begin to understand this interface, we 
allow for a reconstruction by placing a SiC bilayer plus an interface containing up to three 
additional atomic layers between the bulk and the multi-layer graphene film [see Fig. Q]. We 
then write the interface structure factor as: 

5 

Fi(£)=J2f^)P^ eZj/CSiC ^ ( 2 ) 

where pj is the relative atom density for the j th interface layer (pj = 1 for a bulk layer 
corresponding to 8.22 x 10 _16 atoms cm -2 ) at a vertical position zj (the zero height is chosen 
as the last layer of atoms in the interface). fj(£) is the atomic form factor of C or Si. The 
fifth atom layer is added to explore the possibility of adatoms between the SiC and the 
graphene. 

To be completely general the scattered amplitude from the graphene film takes into 
account the possibility of a lateral distribution of varying graphene layers. This is done by 
defining an occupancy parameter p n as the fractional surface area covered by all graphene 
islands that are n graphene layers thick. p n is subject to the constraint equation J^Pn — 1; 
where p$ is the fraction of area that has no graphene. The multilayer graphene structure 
factor can then be written in the general form: 



F G (£) 



ilZm/cQ 



(3a) 



n=l lra=l 



m < 2 



D + D 1 + (m- 2)D G m > 2 



(3b) 
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fc is the atomic form factor for carbon and iV max is the number of layers in the thickest 
graphene film on the surface. D is the spacing between the bottom layer of an island and 
the last atom layer in the interface. D\ is the spacing between graphene layer- 1 and layer-2, 
while Dq is the average layer spacing between graphene in subsequent layers [see Fig. [1] . 
Because STM studies of multi-layer graphene films grown on the Si-face indicate some 



buckling of the graphite layer, 19j, |25|, |26| we must allow for a small vertical height distri- 



bution in each graphene layer that gives rise to a structure factor, F m (£) in Eq. fl3a|) . for 
each layer. A vertical modulation of the graphene layers can be modeled two ways. The 
simplest method is to assume an average, layer independent, random vertical disorder, ctq, 
that will give rise to a Debye- Waller term for each layer, i.e. F m {i) = e -9 ^ "^/ 2 . Because 
the vertical modulation is known to decay quickly after the first graphene layer, [25] a more 
refined model uses the same average Debye- Waller term for the upper graphene layers but 
allows for a different distribution of carbon atoms in the first graphene layer such as those 
calculated in ab intio calculations. 48j In this case the the structure factor of the first layer, 
F\(£), needs to be known. As we'll show, both models give very similar results. 

Reflectivity data for a Si-face multi-layer graphene film is shown in Fig. [2j The main 
bulk 4H-SiC peaks occur at £ = 4 and £ = 8. The sharp peaks at £ = 2, 6 and 10 are the 
"quasi-forbidden" reflections of bulk SiC.[46j In SiC reciprocal lattice units, the graphite 
bulk reflections are nominally expected at £~3, 6 and 9 (i.e., ^ = ^G( c Sic/ c G) ; where £q = 
0, 2, 4 etc.). 

We have tested a number of structural models for the graphene/4H-SiC(0001) interface. 
While the majority of experimental studies point to a complicated interface structure, simple 
models consisting of a nearly bulk terminated substrate with the graphene on top are still 
being proposed. [50] However, such models always give poor fits to the x-ray data. This is 
demonstrated in Fig. [2] where we plot the best fit reflectivity for a bulk terminated surface 
where the 6\/3 interface layer is essentially a graphene layer as suggested by Emtsev et 
al..[s(3] In this model, the interface layer-0 in Fig. [1] is replaced by a single carbon layer 
with a graphitic density. The last bulk bilayer density is kept constant at the bulk value 
while the bilayer spacings are allowed to relax (the relaxation from the bulk value are small; 
A 2 c = — 0.03Aand A 2 si = 0.01 A). All other parameters in Eqs. and ([3]) are allowed to 
vary to achieve the best fit shown in Fig. [2j This includes the distance between the last SiC 
bilayer and the graphitic interface layer, Do, and the distance between this layer and the next 
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FIG. 2: (color online) Specular reflectivity versus gj_ (in r.l.u.) for a graphitized 4H-SiC(0001) 
Si-face surface. Circles are the data. Fits to the two model structures in the figure are given. In 
the schematic models filled and open circles are C and Si atoms, respectively. Shaded circles are 
Si adatoms. Dotted blue line is a fit to a bulk terminated SiC(OOOl) surface with a single relaxed 
bilayer. Solid red line is a fit to a model similar to the relaxed bilayer but with the addition of a 
layer of Si-adatoms with p a d= 0.21. 

graphene layer, D\, that relax to best fit values of 2.55Aand 3.62A), respectively. Note that 
D\ is larger than the bulk graphite spacing of 3.354A. As can be seen in Fig. [2j this model 
give a very poor fit to the data at values of t = 5 and 9, the anti-Bragg points for SiC. This 
is typical of all bulk terminated models including those with a substantial modulation of the 
first graphene layer such as the calculated (6\/3x6v^3)R30 surface of Varchon et al. 48 1 and 



Kim et al. 



49J The calculated reflectivity from this theoretical interface, including a relaxed 



bulk surface and the structure factor Fi (£) of the rippled first graphene layer, gives similarly 
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poor fits to the reflectivity near the anti-Bragg positions. 

Better fits can be obtained by an extended interface where an additional partial layer of 
adatoms is added to the simple relaxed bilayer model [see the schematic models in Fig. [2] . 
As demonstrated in Fig. [2j the additional density from the adatoms begins to correct many 
of the deficiencies in the relaxed bulk model fit at the anti-Bragg points (especially near 
£ = 9). Note that the adatom model used to fit the reflectivity is very similar to the model 



261 ] . including the Si adatom density which is p a( j = 0.21 compared 



proposed by Rutter et al. 
to 0.22 in their model. 

The improvement in the calculated intensities by adding an adatom layer is due to the 
increased scattered intensity at the SiC anti-Bragg condition that would be zero for a bulk 
terminated interface. Regardless, the simple adatom model cannot reproduce a number of 
features in the reflectivity data for £>4. The inability of this model to fit the experimental 
data is a result of both an insufficient atomic density in the interface and the atomic gradient 
through the interface. Increasing the width of the interface adds an additional Fourier 
component in Eq. (J2J) that both broadens the fit near I ~ 6 and removes the interference 
minimum at i ~ 9. Therefore, to improve the fits, it is necessary to change both the 
atom distribution and the thickness of the interface layer. The need for an additional 
plane of atoms is also consistent with number of STM experiments of the SiC(OOOl) (6\/3x 
6V3)R30 interface. STM images of "trimer-like" structures suggest at least one additional 
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34 ] 



partial layer of atoms. 

Adding a 5 th layer of atoms and changing the atom density in the interface layer-0 leads 
to a set of nearly identical structures that are shown in Fig. [3j These structures are distin- 
guished by whether an atom plane in the interface is composed of either carbon or silicon 
atoms. This is because in x-ray diffraction the ratio of atomic form factors of Si and C used 
in Eq. ([2]) is determined, to first order, by the ratio of their atomic numbers 14/6 = 2.33. 
Therefore, the model calculation should give a similar fit if all silicon atoms are replaced 
by carbon atoms with 2.33 times the density (however, there is a substantial difference in 
the I dependence of the Si and C atomic form factors that affects both the p's and z/s in 
the final fits). While x-ray reflectivity data alone is unable to distinguish between different 
silicon and carbon compositions in the interface layer, spectroscopic data from a number of 
experimental groups place restrictions on the atomic makeup of layer-0 in Fig. [Tl 

Angle resolved photoemission spectroscopy (ARPES) studies by Emtsev et al.|5jj as well 



10 



(a) 
Bulk 



p =1.0 

P, a =1-0 
p =1.0 



=1.0 O 



qQ 0.63A 

Q= 



1.89A 
0.63A 



1 Graphene Layer 



1 Graphene Layer 



2.31A 




Interface 232A 
layer 




(b) 



Si-up 
adatom 



p,«,=0.82 <T) 
£=0.86 O 



2 32A ^ Graphene Layer 




■0.89A 



c sas-rj: 

adatom Pu =0Jl ° 

P«=0.81 ^O- 
fe=0.85O— 



(d) 



adatom 



Interface 



r °' 86A J layer 

0.75A 
1.61A 



m Bulky 



FIG. 3: Three graphene/SiC interface ball models for graphene grown on the Si-face of SiC deter- 
mined by surface x-ray reflectivity; (b) Si-up model, (c) Si-down model and (d) C-adatom model. 
Open circles are silicon atoms and shaded circles are carbon atoms. The densities p are relative to 
the densities of bulk SiC shown in (a). 



as x-ray photoemission spectroscopy (XPS) studies by Johansson et al.|23j conclude that 
the interface layer has a significant carbon concentration (at least 1.3 times more than the 



carbon in a bulk SiC bilayer 



231 ]) that rules out a purely silicon interface. In fact the x-ray 



diffraction also rules out a purely silicon interface because the density of silicon required 
to get reasonable fits to the reflectivity data is almost half the density required for an sp 3 
silicon film, which is physically unreasonable. These spectroscopic constraints reduce the 
number of possible structures that are compatible with the reflectivity data to three: the 
"C-adatom" , "Si-up" , and "Si-down" models shown in Fig. [31 
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FIG. 4: (color online) Specular reflectivity versus q±_ (in r.l.u.) for a graphitized 4H-SiC(0001) 
Si-face surface. Circles are the data. Solid black line is the best fit to the Si-up model in Fig. [31(b) 
with Dq = 2.3lA(the C-adatom and Si-down model fits would be indistinguishable). Dashed red 
line is the same fit but with a smaller Do = 2.13A. Dotted blue line is the same fit but with a larger 
D = 2.53A. 

The best reflectivity fit to the data is nearly identical for all three models and is shown 
in Fig. HI Table [I] gives the fitting parameters for all models (uncertainty limits include 
variations from sample to sample). In the C-adatom model a carbon rich layer composed of 
three carbon layers is sandwiched between the graphene and a distorted SiC bilayer. The 
total density of these three interface layers is p = 0.61+1.38+0.71 = 2. 70± 0.15. This density 
is lower than the density of a graphene sheet (pc — 3.13) but is 30% higher than a bulk 
SiC bilayer. The two Si models are similar to the C-adatom model in that they contain a 
carbon rich layer, although it is composed of two rather than three carbon layers sandwiched 
between the graphene. The total carbon density of the interface layers in the Si-up and Si- 
down models are p = 1.44 + 0.74 = 2.18+ 0.15, and 0.58 + 1.50 = 2.08+ 0.15, respectively. 
These values are similar to the total bulk bilayer density (p = 2.0) needed to form a sp 3 
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TABLE I: Best-fit structural parameters for graphite covered 4H-SiC(0001) Si-face. Data for both 
the "Si-up", "Si-down" and "C-adatom" models give nearly identical fits. Parameters are defined 
in Fig. [Tj 





dnrl (A) 


Pad (A) 


Aifc(A) 


Pih (A) 


Ala (A) 


Pla 


A 2 c(A) 




A 2 si(A) 




Si-up 


0.88 


0.23 


-0.23 


1.44 


-0.47 


0.74 


-0.29 


0.86 


-0.16 


0.81 


Atom Type 




silicon 




carbon 




carbon 










Si- down 


0.91 


0.58 


-0.28 


1.50 


-0.48 


0.26 


-0.28 


0.86 


-0.16 


0.82 


Atom Type 




carbon 




carbon 




silicon 










C-adatom 


0.89 


0.61 


-0.25 


1.38 


-0.48 


0.71 


-0.28 


0.85 


-0.15 


0.81 


Atom Type 




carbon 




carbon 




carbon 










uncertainty 


0.04 


0.08 


0.04 


0.08 


0.04 


0.1 


0.05 


0.1 


0.05 


0.08 



bonded carbon layer. The two Si models are distinguished by a low density of Si atoms 
either atop or below the carbon rich interface layer. 

There are two similarities between all three models. First, the high carbon densities in all 
three models suggest a complicated carbon bonding geometry that is neither like bulk SiC 
or like graphene. This has also been noted by both Emtsev et al. 5CJ, |5J| and Johansson et 



al. 



231 ] who studied the (6\/3x6v^3)R30 reconstruction that forms before graphitization and 



is known to persist after a true graphene layer has formed. The ARPES data of Emtsev et 
al. 5o[ IotI | shows that the interface layer has a bands (although shifted to higher energy) but 
no 7r bands. This suggests that the carbon concentration is high enough to at least loc ally 
support a sp 2 bonding geometry. In addition, the XPS studies of both Johansson et al. [23J 



and Emtsev et al. 



511 ] find two surface related C Is core level shifts. As a result both studies 



conclude that the (6a/3x 6a/3)R30 surface contains a large amount of non-graphitic carbon 
in inequivalent surface sites in spite of the a bands. 

The second similarity between these models is that, unlike the SiC(0001) C-face, 4lf] 
the Si-face interface reconstruction extends deeper into the bulk. The bilayer between the 
interface layer and the bulk is substantially altered from a bulk bilayer in both density and 
bonding. This deep reconstruction on the (0001) surface is consistent with the prediction 



of Johansson et al. 



231 ] based on relative intensity ratios of surface to bulk XPS peaks. 



Regardless of the model, the bond lengths between this bilayer and both the bulk and the 
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TABLE II: Structural parameters for graphene grown in UHV on 4H-SiC(0001) Si-face (determined 



from this work) and those from furnace grown 4H-SiC(0001) C-face graphene measured in Ref. 411 ] . 
Parameters are defined in Fig. [T] 





Do (A) 


Di (A) 


A?(A) 


°G (A) 


7SiC 


Si-Face 


2.32 ±0.08 


3.50 ±0.05 


3.35 ±0.01 


0.16(-.05/ + .02) 


0.7 ±0.1 


C-face 


1.62 ±0.08 


3.41 ± 0.04 


3.368 ± 0.005 


< 0.05 


0.03 ±0.01 



interface layer are contracted by ~ 17% from the bulk Si-C bond length making them similar 
to the bonds in diamond (1.54A). 52]. 



IV. DISCUSSION 



The main difference between the "C-adatom" and either the "Si-up" or "Si-down" models 
is the low density Si layer in the interface region. While x-ray data alone cannot discrimi- 
nate between an all carbon interface and a carbon rich interface with silicon, spectroscopic 
measurements strongly favor the two models with silicon in the interface. XPS and pho- 
toemission spectroscopy (PES) experiments conclude that, after graphene has formed, a 



significant fraction of Si remains at the interface 
hansson et al. 



23 



34 



531 ] A complete XPS study by Jo- 



'ace-related C Is core level shifts, there 



541 ] This is consistent with a Si adatom 



231 ] find that in addition to two sur 
are also two surface-related Si 2p core level shifts. I 
layer in the interface and a modified Si-C bond between the interface and the bulk-like bi- 
layer below. We note that the bonding configuration of the Si-up model is very similar to a 
model proposed from STM images of the (6v / 3x6y / 3)R30 interface structure below a layer of 
graphene. [261] The model of Rutter et al. 26J suggests an adatom density of p a d = 0.22; within 
error bars of the x-ray value in Tabled] At the moment, however, there is no experimental 
data that can exclude either of the two Si models. 

It is worth pointing out a number of important structural differences between Si-face 
grown graphene and C-face grown graphene. Table [III shows a comparison of structural 
parameters determined from x-ray reflectivity data for graphene grown on the SiC(0001) 
and (0001) surfaces. The distance between the first graphene layer and the interface for 
Si-face graphene film is D = 2.3 ± 0.08A. Figure H] shows the sensitivity of the fit to either 
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increasing or decreasing the value of D by 9%. The measured value of D is large compared 
to the bilayer distance in bulk SiC (1.89A) and at the same time less than the graphite 
interplanar spacing of 3.3541,0 Note that the best fit value of Dq is similar to the value 
of 2.5A measured by STM.[26[] The Si-face value of Dq is larger than the value for furnace 



grown graphene on the SiC(OOOl) C-face, implying that graphene is more strong 



to the C-face interface. This is consistent with ab initio electronic calculations. 



36 



y bound 



371 ] and 



the conclusions of previous inverse photoemission, PES and XPS studies [13j, [14], [55j (recent 
ARPES experiments have a contradictory interpretation to these studies, suggesting that the 
C-face grap hene grown in UHV is in fact less tightly bound to the interface compared to the 
Si-face 50J). It is also worth pointing out that Table HT1 compares Si- face graphene to C-face 
graphene grown in a furnace; UHV C-face graphene used in previous studies [lJa, 114], |50J grows 



nearly 20CPC lower in UHV compared to graphene grown in a furnace. [17] This temperature 
difference may influence the interface structure or order and should be considered when 
comparing results from different groups. 

In addition to Dq, Table HT1 also shows that the spacing, D%, between between graphene 
layer- 1 and layer-2 is larger than the spacing between all other graphene layers (for both 
Si-face and C-face graphene). Note that the error bar on Dq is significantly larger than 
those reported for C-face graphene films. This is because furnace grown C-face graphene 
films are much thicker than UHV grown Si-face films. Thinner films broaden the graphene 
Bragg peaks at £ = 3, 6, and 9, making the peak position and thus the layer spacing more 
difficult to measure. Table [U also shows that the roughness of the SiC surface is more than 
an order of magnitude larger for Si-face graphene than for C-face graphene (the surface 
roughness is proportional to 7sic)- This is consistent with earlier x-ray results [5j measuring 
growth-induced substrate roughness and with the substrate roughness observed in many 



STM and AFM images. 
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30 



56] 



The most important result of this work is that the interface layer for Si-face graphene 
is not a simple relaxed bulk termination of the SiC surface. This has a bearing on how to 
interpret electronic structure calculations of the graphene-SiC interface. To date, ab initio 
electronic structure calculations of the graphene/SiC(0001) system have started from a flat 



graphene layer placed above an idea 



into a slightly distorted bilayer. 



26 



ized bulk terminated SiC surface that is allowed to relax 



36 



371 ] These calculations use an artificially contracted 



graphene sheet that is commensurate with a small SiC cell to allow for reasonable fast 
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calculation times. Rutter et al. 



261 ] have looked at a Si adatom model but, as with other 



calculations, used a simple relaxed bulk SiC bilayer below the adatoms. The result of all 
these calculations is that the first graphene layer above the relaxed bulk SiC bilayer acts 
as buffer that partially isolates the electronic properties of the next graphene layer from 
the substrate. 36|, |37| While these calculations are an important first step in predicting the 
existence of a buffer layer, their ability to predict the structure of the interface and thus its 
electronic properties is a concern given that the x-ray results show a much more substantial 
reconstruction that has few characteristics of a bulk SiC bilayer. Further experimental 
evidence for a buffer layer comes from ARPES, where a carbon-rich layer with substantial 



,3 



sp 2 -bonding is found without any indication of it bands characteristic of graphene. 23 
ARPES measurements also clearly show 7r-bonded graphene layers above this carbon-rich 
layer, although there are different interpretations of spectral structure very close to the K- 
point of the graphene Brillouin zone for the first of these structural graphene layers. 57|, [58] 



We suggest that in the ab initio results, the first graphene layer in the calculation mimics, 
to some extent, the properties of the interface layer-0. In fact recent ab initio calculations, 
using a full (6\/3 x 6v^3)R30 cell, find that the first graphene layer is significantly distorted 
from a flat graphene sheet and does not show the dispersion characteristic of an isolated 
graphene sheet. 48|, |49J In fact the calculated modulation amplitude is 1.23A, a value not 



far from the ~1.8A interface layer width measured for all the structural models in Fig. [31 
These results suggest that the assumption of a distorted thick carbon-rich layer acting like 
a buffer layer may be correct. A more realistic interface calculation will be necessary to test 
this assertion. 

In addition to the interface structure, the reflectivity data provides additional evidence 
supporting conclusions based on STM and LEED that the (6v^3 x 6v^3)R30 reconstruction 
observed after graphitization is a true reconstruction of the graphene film. Since the discov- 
ery of the (6v^3 x 6v^3)R30 LEED pattern, it has been suggest ed that it is a moire pattern 



due to the near commensuration of graphene with SiC. 10|, |2Cj, 21 1 Early STM experiments 



supported this claim because they imaged a 6x6 reconstruction |19l. 1201. 1211. 1531. |59| instead of 
the (6v / 3x6v / 3)R30 pattern observed in LEED. More recent STM experiments on the other 
hand, have directly imaged the (6\/3 x 6y / 3)R30 structure and shown that the graphene 



has a vertical modulation with this lateral periodicity. 25j, |26( Recent ab initio calculations, 
using a full (6\/3 x 6v / 3)R30 cell, find that at least the first graphene layer above the SiC 
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has a substantial modulation amplitude. 48j The x-ray reflectivity data supports these re- 
cent experiments and confirms that there is a vertical modulation of graphene grown on the 
Si-face. This can be seen by comparing the C-face and Si-face graphene layer roughness or 
corrugation, an, from Eq. ( l3a|) [see Table HI] . oo is much larger on Si-face grown graphene 
than on C-face furnace grown graphene. gq is determined almost solely by the intensity 
decay of the graphite Bragg points as a function of I. Because of the exponential form in 
Eq. fl3aj) . a finite Oa manifests itself as a decay in the graphite Bragg peak intensities at £=6 
and 9. This is demonstrated in Fig. [5] where we compare fits from a graphene film with no 
roughness (<tg = 0A) to a fit with a large roughness (<jg = 0-3A) and to the best fit value for 
the Si-up mode in these experiments <tg = 0.16A. Note in Table HT1 that a n is essentially zero 

nrinn 

for C-face graphene films that show no LEED reconstruction patterns, [5|, HO, H3|, HJj We can 
interpret Oa as originating from an actual modulation of the graphene film, but the value of 
0.16A is considerably smaller than the value of 0.6A measured by STM for the first graphene 
sheet above the interface layer. 26] This difference arises in part because STM is measuring 
a modulation in the electron density of states instead of an actual structural modulation 
but more importantly, because the x-ray value is an average over all the graphene layers 
in the film. Riedl et al.[25| have shown that the vertical modulation amplitude decays by 
approximately a factor of two from the first to the second graphene layer. Therefore, thicker 
graphene films weight the average gq to lower values. 

We can estimate the first graphene layer modulation, Oq , from the measured mean 
modulation <jg if we assume that the modulation decays in subsequent layers as = 
<7^q exp|— \D n] w here XD = ln(2) (the factor of ln(2) assumes the amplitude decay measured 
by Riedl et al. 25] is correct). To calculate we only need to know the relative amount of 
graphene that is thicker than N layers, P^. P/v is calculated from the areal coverage, p n 's 
in Eq. (I3a| ; P/v = C Yln^N Pn (C is a normalization constant). Then is given by; 

= °G P nl Yl P " [~ XDn \ ■ ( 4 ) 

n n 

The measured distribution of graphene layer thickness, p n , for a nominally 2-layer graphene 
film is shown in the layer height histogram in Fig. HVl The average number of graphene 
layers is 1.9 ±1.5. The distribution is very wide, in part reflecting the spatial average 
over the large x-ray beam footprint (the footprint is bigger than the sample width of 3mm 
when £< 1.8). In particular the high areal fraction not covered by graphene (18%) can be 
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associated with slow growth kinetics at cooler substrate regions near the edge of the sample 
caused by non-uniformity in the e-beam heating. These wide distributions are also seen in 
low energy electron microscope (LEEM) images. 3j], 32] Using this measured distribution, 



we estimate to range between 0.5-0.8Areflecting the uncertainty in the measured value 
of og- This result is in very good agreement with the STM value. 26 1 

For comparison, we can estimate the vertical modulation of the first graphene layer using 
a different method. This is done by using a modulated first graphene layer on-top of the 
Si-up adatom model to calculate F\{£) in Eq. (l3al) . To do this we use the graphene structural 
coordinates calculated by Varchon et al. [481] The modulation of higher graphene layers are 
included using the rms roughness in the standard model. While the calculated relative 
vertical carbon positions are maintained in the first graphene layer, the absolute positions 
are scaled by a multiplicative constant so that the peak-to-peak amplitude can be varied. 
The best fit structure to the reflectivity data gives the first graphene layer amplitude to be 
0.82Apeak-to-peak, which is slightly larger than than the range of (Jq estimated above. 



It is worth comparing the 



the simpler AES method. 11 



grap 



60, 



rene thickness measured by SXRD and an estimate from 



61] The ratio of the Si(LVV)/C(KLL) peak area can be 



used, along with proper electron mean free paths and excitation cross sections, to esti- 
mate of the number of graphene layers provided a suitable model for the interface is used. 
Most grou ps use an interface model consisting of a bulk terminated substrate with graphene 



above. 



11 



601 ] (Ref. [6l| considers a \/3 x Si adlayer). While bulk termination is clearly 
inconsistent with both models in Fig. [31 it can still be used with the proviso that it will 
overestimate the film thickness by ~ 1 layer because the measured C(KLL) intensity includes 
a contribution from the dense non-graphitic interface carbon layer in Fig. [31 AES measure- 
ments on the same sample as the data for Fig. IIVI estimate the average graphene coverage to 
be 3.2 layers compared to 1.9 by x-rays, consistent with the lack of a realistic interface layer 
in the AES calculation. Because the electron penetration depth of the Si(LVV) electron is 
short, AES estimates become more uncertain for graphene layers exceeding four layers. This 
makes the AES method more applicable to Si-face films than to thick furnace grown C-face 
graphene films. 
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FIG. 5: (color online) Specular reflectivity versus qj_ (in r.l.u.) for a graphitized 4H-SiC(0001) 
Si-face surface. Circles are the data. Solid black line is the best fit to the Si-up model with 
o"g = 0.16A. Dashed red line is the same fit but with o"g = 0.30A. Dotted blue line is the same fit 
but with a larger <tg = O.oA. 

V. CONCLUSION 



We have used surface x-ray reflectivity measurements to determine the graphene/SiC 
interface structure for graphene grown on the 4H-SiC(0001) Si terminated surface (Si-face). 
We find that the interface is not composed of a simple graphene-like layer above a relaxed 
SiC bilayer that has been recently proposed. [5^] Instead, the interface reconstruction is 
more complicated and extends deeper into the bulk. Based on these x-ray experiments and 
previous XPS studies 2iJ the best model for the interface is composed of a substantially 
relaxed SiC bilayer, above which, a dense carbon layer containing a partial layer of Si 
atoms separates it from the graphene film. The carbon density in this intermediate layer is 



approximate 
STM 
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25. 



y 2.1 times larger than in a SiC bilayer. This model is consistent with previous 



experiments. 



261 ] and XPS 23( results as well as the lack of interface ir bands in ARPES 



50. 



5l| The bond distance between the Si adatom layer and the first graphene 
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FIG. 6: The normalized probability, p n , of a n-graphene layer stack from a UHV grown Si- face film 
as determined by x-ray reflectivity. The x-ray average is 1.9±1.5 graphene layers while the AES 
estimate of the average is approximately one layer thicker (3.2 layers). 



layer (Si-up model) is 2.32±0.08A. While this distance is short compared to the interplanar 
graphene spacing, it is still larger than the corresponding distance measured on furnace 
grown SiC(OOOl) C-face graphene (1.62 A), 41] indicating that the graphene on the Si-face 
is less tightly bound to the substrate than furnace grown C-face graphene. This difference 
is consistent with ab-initio calculations. 36J] We propose that the dense carbon layer with Si 



adatoms plays the role of the buffer layer predicted by Varchon et al.|36] It is this layer that 
partly isolates subsequent graphene layers from interactions with the substrate. 

The x-ray results also indicate that the graphene is not flat but has a corrugation 



amplitude consistent with that measured by STM.[25|, |26| This result confirms that the 
(6V3x6V3)R30 periodicity imaged by STM is dominated by a real structural corrugation 
of the graphene and not by changes in the local density of states. 
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